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Aristotelian electrodynamics (AE) describes the regime of a plasma with a very strong electric field that is 
not shorted out, with charge current determined completely hy pair production and the balance of Lorentz 4- 
force against curvature radiation reaction. Here it is shown how the principal null directions and associated 
eigenvalues of the field tensor govern AE, and how force-free electrodynamics arises smoothly from AE when 
the eigenvalues (and therefore the electric field in some frame) vanish. A criterion for validity of AE is proposed 
in terms of a pair of “field curvature scalars” formed from the first derivative of the principal null directions. 


I. INTRODUCTION 

Gruzinov has proposed that force-free electrodynamics is 
inadequate for describing some relativistic magnetospheres, 
because of insufficient pair supply to short out the induced 
electric fields. In particular he argued that this is the case for 
“weak pulsars” [1]. In this case, parts of the magnetosphere 
would fail to be force-free. When the fields are very strong 
and E • B is nonzero, charges will be accelerated to near the 
speed of light, reaching a steady state where the work done 
on them by the electric field is transferred to their radiation 
fields. In this situation, all charges move approximately at 
the speed of light, and the current is completely determined 
by the supply of charge. This regime was called Aristotelian 
electrodynamics (AE), because it is the velocity rather than 
the acceleration of the charges that is determined by the field. 
The name is quite fitting, since a central focus of Aristotelian 
mechanics was the regime of terminal velocity, in which the 
transients can be neglected [2]. Consideration of this regime 
in pulsar electrodynamics was earlier discussed in [3, 4]. 

In the degenerate case when E • B = 0, if the field is mag¬ 
netic {B^ > E^), there exist Lorentz frames with vanishing 
electric field. The charges therefore have no secular accel¬ 
eration, so it can be a good approximation in strong, mag¬ 
netic degenerate fields to neglect the energy and momentum 
of the charges altogether. In this, force-free approximation, the 
current is entirely determined by the fields at one time (see 
e.g. [5] for a review and references). The Aristotelian and 
force-free approximations may coexist in one physical sys¬ 
tem that includes both degenerate and non-degenerate regions. 
Both approximations result from neglecting of the mass of the 
charges in certain quantities. Together they comprise what 
Gruzinov called the electrodynamics of massless charges [6]. 

Besides pulsars, another astrophysical setting where one 
would expect AE conditions to prevail is the magnetosphere 
of a compact binary system near coalescence. Even with co¬ 
pious pair creation, because of the motion of the binary, the 
fields change too rapidly for the charges to short out the elec¬ 
tric field. One would therefore expect very strong charge ac¬ 
celeration to take place [7]. AE may also be relevant in a 
laboratory setting, for plasmas in Petawatt laser fields [8]. 
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In this paper I reformulate AE in a relativistically covari¬ 
ant form, and aim to elucidate the structure of the theory and 
its interface with force-free electrodynamics (EEE). The key 
observation is the central role played by the principal null 
directions of the electromagnetic field. Beyond that there is 
nothing new here. The covariant form also applies in curved 
spacetime, so would be particularly useful in systems contain¬ 
ing a black hole. The spacetime signature is (H-) and the 

speed of light is sometimes set to c = 1. 


II. STRUCTURE OF ELECTROMAGNETIC FIELDS 

A non-null electromagnetic field Fab at a point has two null 
eigenvectors k±, with opposite eigenvalues, 

F\k^^=±Eokl. (1) 

The directions of these null eigenvectors are called the prin¬ 
cipal null directions (PNDs) of the field [9]. The eigenvalue 
EEq is the electric field in a frame in which the electric and 
magnetic fields are parallel or one of them vanishes. In AE, 
Eq f 0, and positive charges travel along one PND while neg¬ 
ative charges travel along the other. In PEE, Eq = 0, and the 
total current 4-vector also lies in the timelike plane spanned 
by the two PNDs. Eor null fields there is only one null eigen¬ 
vector, with vanishing eigenvalue. 

A simple way to understand this structure is to observe that, 
at a given spacetime point, any electromagnetic field 2-form 
F can be presented in one of two canonical ways in terms of 
an adapted orthonormal set of 1-fotms {dt, dx, dy, dz}: 

^generic = Eq dz A dt + Bq dx A dy, (2) 

= FQdzA{dt-dx). (3) 

The field in (2) has two null eigenvectors, k± = dt Ed z, with 
eigenvalues EEq. The field in (3) has one null eigenvector, 
dt + dx, with vanishing eigenvalue. The electric and magnetic 
fields in (2) are parallel in the Lorentz frame defined by dt- 
The general statement is that they are parallel in any frame 
lying in the timelike plane spanned by fc_(_ and k_. If the field 
is degenerate, i.e. if F A F = 0 (E • B = 0), then either 
Eq = 0 or i?o = 0; that is, either E or B will vanish in 
these “field eigenframes”. If the field is null then E and B are 
perpendicular and have the same magnitude in all frames. 
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FIG. 1. The space of electromagnetic fields modulo Lorentz trans¬ 
formations. The coordinates {Bq, Eq) are the invariants defined in 
(4,5). Degenerate fields are represented by the two axes. (±f3o, 0) 
are identified on the dashed lower boundary, so the space is actually 
a cone. AE applies above the grey strip of height where R 

is the curvature radius of the particle paths, reckoned in an instan¬ 
taneous field eigenframe, and units e = me? = 1 are used. The 
dashed Bq axis corresponds to magnetic or null {Bq = 0) degener¬ 
ate fields; FFE applies for these if the energy density of the field is 
much greater than that of the charges. 


The quantities Eq and Bq are related to Lorentz scalars by 

- Bl = E^-B^ = -^EabF‘^^ (4) 

EqBq = E • B = (5) 

To fix signs uniquely, we take Eq > 0. Eq and Bq are given in 
terms of a = E^ — B^ and b = 2E • B by (y/+ 6^ ± a)/2, 
respectively. The PNDs can be specified by a pair of spatial 
unit vectors v± in a given frame via k'^ -o- (1, v±). Gruzinov 
gives an explicit expression for v± in terms of the electric and 
magnetic fields: 


v± = 


E X B ± {EqEi + BqB) 


F? 


52 


( 6 ) 


[The denominator can also be written more symmetrically as 

1(52 + 52+52 + 52).] 

The PNDs can be constructed as quadratic expressions from 
the spinor factors (or eigenspinors) of the electromagnetic 
spinor [9].' The Newman-Penrose formalism or other spinor 
techniques might therefore be useful for analytical and/or nu¬ 
merical studies of AE. 

The set of fields modulo Lorentz transformations can be 
represented as the upper half plane with vertical axis Eq and 
horizontal axis Bq (see Fig. 1). Degenerate fields have either 
5o = 0 or 5o = 0, while null fields have Eq = Bq = 0. 
When Eq = 0, the two values ±5o label the same set of 


fields. The moduli space of fields is thus the cone obtained 
by identifying the positive-5o half of the dashed line in Fig. 1 
with the negative half; the null fields lie at the vertex, and the 
magnetic and electric degenerate fields correspond to a pair of 
opposite rays on the cone. 

The current 4-vector lies in the plane spanned by the PNDs 
in both the AE and FFE regimes, though it is determined by 
different conditions. The distinguishing factor between these 
regimes is whether the PND eigenvalue is nonzero or zero. In 
the following we characterize these two regimes. 

The PND’s define a pair of “field curvature scalars”, 


,-l _ |(fc± ■ V)fc=Fl 
" \k+-k_\/2 ■ 


(7) 


These quantities have dimensions of inverse length and are 
invariant under arbitrary rescalings of k± ? The two radii R± 
coincide if the PNDs are surface forming. For static electric 
or magnetic fields they coincide and are equal to the curva¬ 
ture radius of the field lines. More generally, they measure a 
rate at which one eigendirection bends away in the orthogonal 
spacelike direction when moving along the other PND. They 
are equal to twice the magnitudes of the spin coefficients tt 
and r associated with a null tetrad constructed from the PNDs 
[9]. They may serve generally to define the relevant length 
scales determining the applicability of the AE approximation 
and the “guiding center approximation” (see e.g. [10] and ref¬ 
erences therein).^ 


III. FORCE-FREE ELECTRODYNAMICS 

If Eq = 0, then the field is magnetic or null degenerate. If 
the energy density of the field is much greater than that of the 
charges, then the energy and momentum of the charges can be 
neglected, which means that the 4-force on the current can be 
set to zero, 


Fabf = 0 . ( 8 ) 

This condition in turn implies that the field is degenerate. In 
the magnetic case, the FF current must be a linear combina¬ 
tion of the vectors k^, since the latter span the kernel (i.e. the 
null eigenspace) of Fab- In the null case the PNDs coincide, 
the kernel of Fat is spanned by the unique PND and an or¬ 
thogonal spacelike direction, and the current must be a linear 
combination of vectors in those directions. If there is only one 
sign of charge present, the current cannot be spacelike, so in 
the case of null fields it must be null. 


^ An electromagnetic field can be represented by the trace-free 2x2 matrix 
$ = (E -h iB) • cr. The Lorentz group acts by similarity transformation 
on so the eigenvalues of $ are Lorentz invariant. The properties of the 
Pauli matrices imply = (E + zB) ■ (E -h zB) = [Eq iBo)^I, so 
the eigenvalues are di(-Eo -|- iBo). The eigenspinors X± are the square 
roots of the PNDs of Fab, the sense that o (1, A^crA±), with the 
normalization A^ A= 1. 


^ Since = 0, (A:+ ■ V)/c_ is orthogonal to k—, so (/c+ • V)/c_ = 
Ak— -|- V for some spacelike vector v orthogonal to both k- and /c+, and 
|(/c+ • V)/c_| = |z;|. Under a rescaling k± a±k± the coefficient A 
is modified, and v olj^ol—v, so |z;| is simply multiplied by |a+a_|, 
which cancels against a similar factor from the denominator of (7). 

^ I am grateful to Antony Speranza for suggesting that the invariants (7) 
might supply the relevant notion of curvature length scale in this context. 
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In either case, somewhat surprisingly, the current is deter¬ 
mined by the fields at one time when ( 8 ) and Maxwell’s equa¬ 
tions are imposed. One therefore has stand-alone evolution 
equations for magnetically dominated or null force-free fields, 
without reference to the charge or current densities. In the 
magnetic case the initial value problem for these equations is 
well-posed. [11-13]. 

IV. ARISTOTELIAN ELECTRODYNAMICS 

For any field with Eq ^ 0, i.e. which is not magnetic or 
null degenerate, the electric field can accelerate charges with¬ 
out deflection along the spatial directions of the PNDs. In a 
strong field the charges accelerate until they radiate energy 
at the same rate as the Lorentz force supplies it to them— 
unless collisions with something slows them down before that 
happens.^ The power input is ^ cE^c, reckoned in a field 
eigenframe in which the charge has speed « c, while a charge 
moving with large Lorentz 7 factor on a path following a tra¬ 
jectory with radius of curvature R emits curvature radiation 
with power [15]. Equating the input and output 

power yields 7 = (3£'o^^/2e)^/^.^ The 4-velocity of the 
charge is nearly parallel to the principal null vector if the 
charge is positive, and to if the charge is negative. 

When does AE apply? The electric field should be strong 
enough that it is a good approximation to treat the current of 
the charges as running parallel to the principal null directions. 
If the field is approximately static, it is presumably necessary 
that the 4-velocity is “nearly lightlike” as reckoned in the ap¬ 
proximately static frame, and that the conditions for reach¬ 
ing terminal velocity have been met. Let us use units with 
e = mc^ = 1 for a moment. The conditions just stated 
are that the terminal gamma factor must be relativistic, say 
{EqR^Y^'^ > 2 , and that the particles reach that terminal 
velocity. This latter condition requires that the voltage drop 
over a distance ^ R, reckoned in an instantaneous field eigen¬ 
frame, be greater than the terminal gamma factor times the rest 
energy of the electron, EqR > {EqR?Y/'^, i.e. EqR?/^ > 1 . 
The former condition is then automatically met in any relevant 
situation, since then i? ^ iie = 1 (the classical electron ra¬ 
dius e^/mc^ in these units). This reasoning suggests that the 
condition of AE applicability is 

Eq > (9) 

Note that the range of applicability grows as R grows (see 
Eig. 1). We conjecture that the field curvature scalars (7) pro¬ 
vide the correct general interpretation of the length scale R in 

(9). 

Given a unit timelike vector it“, we can normalize the prin¬ 
cipal null vectors by the condition k'^Ua = 1. Then, because 


'* The charge acceleration could be halted, for example, by inverse Compton 
collisions with ambient photons, as in the stagnation surface scenario for 
black hole jets discussed in [14], 

^ The characteristic frequency of the radiation is ~ 'y^c/R ~ 
(Eo/e)3/‘^ci?-l/2. 


the charge 4-velocities are parallel to fcj in the AE regime, the 
current in a region where Eq ^ 0 takes the form 

r = p+ki-p_k-_, ( 10 ) 

where and — p_ are the densities of positive and negative 
charge in the frame u°‘. In terms of the velocities ( 6 ) relative 
to a fixed Lorentz frame, the 4-current density is given by 

(pJ) = (P+-P-,P-bV+-p-v_). (11) 

Using the eigenvector property (1), the Lorentz 4-force den¬ 
sity acting on the current ( 10 ) is 

= Eo{p+kl + p_k^_). ( 12 ) 

The power deposited into the charges in the frame is 

Eo{p+ + P-). 

The charge density is determined by Eab without time 
derivatives via Gauss’ law, p = V • E. Since the current 4- 
vector ( 10 ) lies in the plane spanned by the two null eigen¬ 
vectors, the remaining freedom in is only one function. If 
all the charges in a given region have the same sign, then the 
current is null and thus fully determined. If instead both signs 
of charge are present (nonzero pair multiplicity), the charge 
densities are determined by pair production and subsequent 
propagation subject to the continuity equations, 

Va(p+fc“) = Va(p_fc“) = r. (13) 

The pair creation rate L depends on Eq and the photon den¬ 
sity. Once the form of L is specified. Maxwell’s equations, 
together with Eqs. (10) and (13) and initial values for p±, de¬ 
termine the time derivatives in terms of the field and charge 
densities at a given time. This system of equations is thus 
naively deterministic. (Whether it defines a well-posed initial 
value problem has not been examined, although Gruzinov has 
evolved them numerically and the solutions seem to behave 
well [ 1 , 16].) 

A. Example: Gruzinov’s device 

To illustrate the workings of AE/EEE, Gruzinov considered 
an arrangement in two spatial dimensions where EEE and AE 
regimes coexist side by side. He called this the “Device” [17] 
(see Eig. 2 in this reference for an illustration). In the Device, 
the y = 0 plane is a conductor, and opposing EEE Poynt- 
ing fluxes propagating in the Ex directions collide in an AE 
zone where the energy is converted to curvature radiation. The 
charges are electrically attracted toward the image charges and 
repelled by the image currents below the conducting plane. In 
the EEE zone the charges move rectilinearly at the speed of 
light, so the electric and magnetic forces balance. This zone 
transitions to an AE radiation zone, in analogy with the tran¬ 
sition outside a weak pulsar, without any discontinuity. In the 
AE zone, the magnetic field is weaker than the electric one, 
so the net force attracts the charges to the conductor. Here I 
use the stationary Device to illustrate the formulation given 
above. 
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The stationary field can be expressed as 

F = E{y)dy f\[dt-\-I3{x)dx\, (14) 

corresponding to an electric field E in the y direction and a 
magnetic field —j3E in the z direction. F is a simple 2-form 
(as is any electromagnetic field in 2 h-1 dimensions), so it is 
degenerate. That is, one or both of Eq and Bq vanish. 

For /3^ < 1 the field is electric, so it is Bq that vanishes. 
There are two PNDs, 

k± = dt- /3dx±y^l- jS^dy, (15) 

with eigenvalues 

Eo = ±E,/1-PF (16) 

The frames in which the magnetic field vanishes are those 
in the plane. From (15) we see that these have x- 

coordinate velocity —/3 and any y-coordinate velocity with 
magnitude less than or equal to -^/l — j3^. A positive charge 
moving in the direction of the null vector k+ feels a Lorentz 
force proportional to fc+. If /3 = 0 the spatial direction of 
such motion is just that of the electric field, dy, while if /3 is 
nonzero there is also a dx component. 

For /3^ = 1 the field is null, and the two PNDs coalesce 
into a single PND with vanishing eigenvalue. In this case the 
field takes the form F = E{y)dy A {dt ± dx), which is a 
solution of the force-free equations. In fact it has the form of 
the simple class of null Poynting flux solutions discussed in 
[5]. For /3^ > 1 the field is magnetic, but it is not a force-free 
solution for any nontrivial choice of the functions E and P, so 
that case plays no role here. 

Let us now consider the AE field equations in the electric 
case /3^ < 1. First, the field satisfies Faraday’s law, dE = 0, 
by inspection. The remaining AE equations require that (i) the 
current computed from F according to Maxwell’s equations 
is equal to (10) and (ii) the continuity equation (13) holds. 
Assuming there is no pair creation (F = 0), and that only 
electrons are present, the continuity equation follows from the 
Maxwell equation, = —p-kt- This equation implies 

/ 9 _ = — E^y, and E^y/E = —/3,x/\/l — Since F depends 
only on y, and depends only on x, the general solution is 
given by 

E = F(0) exp(-y/f), j3 = sin(x/£), (17) 

where £ is a constant length. This solution matches smoothly 
to the force-free Poynting flux solution with /3 = ±1 at 


x|^ = F.'k 12. The charges and Poynting flux are thus flowing 
inwards toward x = 0 from both sides. 

What about the validity of the approximations for the De¬ 
vice? In the FF zone the charges move in straight lines at 
the speed of light, which can presumably be understood as 
an ultrarelativistic approximation. Where the radiation zone 
begins, Eq starts out at zero (16). According to (9), AE 
therefore applies at the transition only if i? = oo, where 
R is the radius of curvature in an instantaneous field eigen- 
frame. We proposed that the invariant i?_ defined in (7) may 
capture the relevant curvature radius. The vector fields (15) 
are surface forming, so = R-, and we find R'^^ = 
9a;(l — This diverges at x/f = 7 r/ 2 , so R± 0 

there, hence for any F( 0 ) there is a region close to x/^ = 
7 r /2 where the applicability of AE is questionable. Specif¬ 
ically, using (16), the condition (9) becomes E(d)e~'^/^ > 
.(“^/^[sin(x/£)/ cos^(x/f)]^/^. 


V. SUMMARY 

Aristotelian Electrodynamics [1] is the theory of a plasma 
with very strong electric field, whose charges move ultrarela- 
tivistically, essentially at the speed of light in the direction of 
the local electric field. In this paper we have formulated AE 
in a way that brings out the central role played by the princi¬ 
pal null directions of the electromagnetic field tensor ( 1 ). 
These determine the current (10) up to the densities of positive 
and negative charge, and can be used to construct the field cur¬ 
vature scalars (7) relevant to applicability of the AE approx¬ 
imation. The total charge density is determined by the diver¬ 
gence of the electric field, so just one function in the 4-current, 
say the positive charge density, remains undetermined by the 
field. This is determined by the pair production rate F and the 
continuity equation (13). The rate of 4-momentum density de¬ 
posited into the charges is proportional to the eigenvalue Eq 
( 1 ), and this eigenvalue can also be used to parametrize F. 
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